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Overview, Objectives and Progress 


Executive Summary 

The Space Engineering Research Center (SERC) at M.LT., started in July 
1988, has completed three and one-half years of research. This Semi-Annual Report 
presents annotated viewgraph material presented at the January 1992 Steering 
Committee and Technical Representative Review. Nine faculty and approximately 
35 graduate students and like number of undergraduates through the 
Undergraduate Research Opportunity Program (UROP) participate in the activities 
of SERC, drawn from the Departments of Aeronautics and Astronautics, Electrical 
Engineering and Computer Science, and Mechanical Engineering. 

The objective of the Space Engineering Research Center is to develop and 
disseminate a unified technology of controlled structures, to codify and disseminate 
this technology, and to train a generation of skilled engineers. The disciplinary 
research in CST continues and in this past year several notable milestones have 
been reached. The interferometer testbed became operational and is an active 
vehicle for CST researchers. The engineering model of the MACE multibody test 
article also was commissioned. It is a pathfinder for control of articulated flexible 
structures. The MODE experiment flew on STS-48. Its goals were to examine the 
dynamics of space structures in zero-gravity. 

Moreover, there has been an ongoing effort toward developmental work in 
CST by the Space Engineering Research Center as we look ahead to the 
university/government/indiistry teaming of technologies for space applications with 
SERC as a model for university/government/industry cooperation. 
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Reference Mission & Science Requirements 
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Schematic of a two dimensional interferometer 
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Reference Spacecraft Design Summary 
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Decision made to address only a subset of the DPL problem. 
Errors due only to internal flexible motion will be considered, that 
due to external flexible and rigid body motion are not measured. 
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The performance metric as defined by internal flexible motion is 
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Disturbance Source & Signal 
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Disturbance Signal (con't): 

- Signal representing RWA spike envelope which increases with 
the square of the frequency up to wheels speed limit at 70 Hz. 
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"MOCK SIDEROSTATS" CAN BE SEEN AS THE THREE PLATES IN THE FOREGROUND. THESE 
E LOCATIONS A, B, AND C AND CONTAIN OPTICAL RETRO-REFLECTOR MOUNTS FOR THE 
SER MEASUREMENT SYSTEM. 
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TWO ACTIVE STRUTS CAN ALSO BE SEEN AT LOCATIONS IN THE TRUSS WORK. 



Interferometer Testbed 
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WE USE A SUN WORKSTATION AS A HOST FOR A REAL-TIME CONTROL COMPUTER 
WHICH IS USED FOR COMPUTING THE PERFORMANCE METRIC AND IMPLEMENTING 
CONTROLLERS. IT INTERFACES WITH THE SIX HP LASER CHANNELS, 16 A/D'S AND 10 D/A'S, 
AND 16 ANTI-ALIASING BESSEL FILTERS. 
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A Tektronix spectrum analyzer is used to drive the disturbance source, to 

MEASURE RELEVANT TRANSFER FUNCTIONS, AND TO HELP MONITOR THE PERFORMANCE 
METRIC. 



Testbed Architecture 
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Passive Damping in the MIT SERC 
Controlled Structures Testbed 
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Viscoelastic Damping 
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Results of Viscoelastic Damping 

Fifty J-struts were added to the structure, placing then to maximize the minimum level of dar 




1 | 

i 




1IHI! R I l!ll Mill! I HI Iffll I III Hill I I ll VMM II I uni nnnniiirimii 1 1 n 1 1 iimi mum ■ r i m >iirin 


Results of Viscoelastic Damping 
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Viscous “Target” Damping 
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D-Strut Model and Component Tests 

The D-strut can be effectively modelled as a lead-lag network. Typical predicted results of the loss factor and 
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D-Strut Model and Component Tests 
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Transfer Function Results of Target Damping 

Consistent with the goal of placing the D-strut viscous dampers for a specific control task, the results ar< 




Transfer Function Results of Target Damping 
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Pathlength Control Using Isolation Mounts 








Pathlength Control Using Isolation Mounts 
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e key results are that J-struts, while offering only a low peak loss factor, have a unit mass and stiffness 







Summary 
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Finite Element Model and Identification 

Procedure 


4> 

CO 

dd 

V 

-4-J 

a; 


kl 

G 

b0 


d 

4> 

d 

CO 

dd 

CO 


d 

t 

u 

4) 

CO 

"D 

"3 

O 

CO 

6 

(m 


4) 

Cu 

6 

d 

-*-> 

d 

4) 

E 


"3 

-4-» 

d 

4> 

+j 

ki 



'2 

CO 

cd 

Q 

d 

4> 


dd 

13 


> 

<-•-1 


o 

o 

-♦-> 

d 


u 

4) 

6 

CO 

u 

U 

d 

J2 

1-4 

a 

CO 

4> 

4 

> 

dd 



’ CO 
CO 

bO 

G 

d 

a 

* V 

u 

d 

O 

"X3 

V 


> 

T3 

s 

» F-1 

-4-J 

U 


V 

dd 


v 

cu 

U 4 

o 

s 

V 

dd 


O 


d 

§ 

cu 

a 

d 

3 


bd 

4> 

>* 

d 

d 

B 

dd 

V 

u 

cu 

<u 

dd 

-4J 

CO 

CO 

d 

u 

co 

"5 

o 


"3 


4) 

co 

O 

cx 

kl 

d 

a 

v 

dd 

H 

-6 

v 

cu 

o 

T3 

> 

4 > 

~o 

v 

dd 


d 

dd 

drf 

3 

CO 

d 

o 

’> 

V 

ki 

CU 


G 

4) 


8 

bO 

d 


4 

Xi 

O 

.2 

4 > 

d2 

CO 

V 

-k> 

H 

C/2 

a 

0 

01 

w 

co 

4) 

dd 


0 

CO 

a 

1 

V 

(H 

d 

a 

a; 

£ 


« 

-4-» 

XJ 

o 

JS 

CO 

d 

d 

d 

'a* 

6 

3 

o 

CO 

-k» 

"O 

o 

JL> 

3 

3 

8 

6 

d 

ki 

V 

4) 

T3 

d3 

d= 

4 

k* 

bO 

o 

d 

d 

-+J 

CO 

rA 

w M 

kl 


v 

JC 


u 

4 


4 

d2 


O 

Oh 

H 

CO 


PRECEDING PAGE BLANK NOT FILMED 




mm i mu min m 




Space 

‘Engineering 

^search 



s 

.2 

03 

u 

ca 

• ^ 

s 

0 ) 

T 3 

t— i 

TJ 

fl 

cd 

13 

TJ 

0 


s 

0 ) 


<D 

H 


0 > 

Sh 

3 

Td 

0 » 0 ) 

•t 3 £ 

d 0 

.a 


a 

r^> 

a 

o 

CO 

5-1 

0 ) 

d 

C 

CD 

S-t 

w 

o 

o 

£ 

cj 

d 

s 

>» 

5-i 

d 

o 

fN 

& 

O to 

S D 

= 

(U 13 
a « 

cc3 <u 


-u 
a 3 

4 w 


d 

d 


bjO 


( 1 ) 

CD 


CD 

0 ) 


CO 

O 

Ph 

H 

CO 


PRECEDING PAGE BLANX NOT FILMED 


22 January, 1992 


Finite Element Model Update 
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Modal Frequency and Damping Comparison 

The plot compares the frequency and damping values for the modes below 100 Hz for the six urouos of transl 
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Procedure effectively fits hundreds of transfer functions. 

Results good enough for control designs. 
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Current Status 



flexibility, inclusion in the FEM lead to improved 
frequency agreement (< 4 % error in first 9) 

fourth vertex stiffened to improve optical alignment, and 
better agreement indicates prior presence of local modes. 
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APPROACH TO THIS PROBLEM RELIES ON THE INTENSIVE USE OF NUMERIC 
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any linear structure can be described by the top equation, a more useful for is 

THE MODAL STATE SPACE FORM SHOWN BELOW. THE STATE VECTOR BECOMES MODAL 
DISPLACEMENT AND FREQUENCY NORMALIZED MODAL VELOCITIES INSTEAD OF PHYSICAL 
DISPLACEMENTS AND VELOCITIES. A USEFUL FEATURE OF THE MODAL FORM IS THAT EACH 
OF THE MATRICES CORRESPONDS EXACTLY TO A DIFFERENT MODAL PROPERTY. 
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USING THESE MATRICES AS A GUIDE, ONE CAN ASSIGN MATHEMATICAL DEFINITIONS TO 
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Single Mass Typical Section 
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THIS EQUATION CAN BE SOLVED IN CLOSED FORM FOR THIS SIMPLE SYSTEM TO YIELD THE 
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Closed Form Solution 
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Asymptotic Behavior of Cost 
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Differences Between Control Levels 
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Closed Form Solution 
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ALREADY USED ON THE SINGLE MASS TYPICAL SECTION. THE NOMINAL DESIGN FOR THIS 
SYSTEM IS A UNIFORM THICKNESS, (1 CM) WITH THE LUMPED MASSES SET TO ZERO. 



Effects of Several Modes 
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Effects of Several Modes (cont) 
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- Design variables: Element thicknesses and lumped mass 
values 

- Mass constraint: Mass must be less than or equal to that of 
a uniform beam of 1cm thickness with no lumped mass 
(nominal) 
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Normalized Subgradients vs. Control 



More actuators than modes reduces importance of control influence 
parameter, p, at high levels of control 
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optimal design if performance is insensitive to disturbances 

In systems where there are more modes than actuators, the effectiveness 
of control at higher control levels will be diminished. 
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A Stochastic Approach to Robust Broadband 
Structural Control 
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Analagous to reflection coefficient 
Reflected power is ( H*H)w*w 
Dissipated power is (/ — 
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Impedance Matching 
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We want to minimize the best estimate of the global performance metric, given that we only have local 
'wledge and local control of the uncertain structure. There are several properties of the structure we can take 
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Control of Stochastic Systems 

These assumptions, and the dereverberated mobility model can now be used to generate a cost functional 
expressed m terms of only the local information, that represents the best guess to the global performance metric. 
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Non-causal impedance match: (maximum dissipation) 


Compensator Design 

First, choose the input power n*,(w) based on the disturbance spectrum V(w) and the dereverberated input 
lhty at the disturbance location G d {u>). Then choose the weighting function C(u>) to approximate the modal 



Compensator Design 
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Choose desired compensator order and use numerical 
optimization. 
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STA Objectives 



MODE: Structural Test Objectives 
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accelerometer locations were chosen so that windowed modes would be 
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Sensor and actuator locations for Straight and L configurations. 




Ground (1 Hz Suspension) vs. Space Torsion Mode (baseline, lowest preload) 
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Ground (1 Hz Suspension) vs. Space Torsion Mode (alpha tight) 
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The last view graph of data was separated into 0-GRAVITY data and an over- 
plot OF 0-GRAVITY AND 1-GRAVITY DATA. THE DATA SHOWN IS FOR THE ALPHA JOINT 
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PRECEDING PAGE BLANK NOT FILMED 


Space 'Engineering %$$tarch Center 


Program Logic 



o 


$ _ 

s 


w 

3 


THE DATA BASE NOW AVAILABLE. 
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This presentation discusses various aspects of the Middeck active Control 
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THE MODE FAMILY OF EXPERIMENTS 



MACE is part of a logical sequence of cost-effective flight 
experiments designed to advance technology of interest to NASA 
in the area of controlled structures. 
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MACE SCIENCE OBJECTIVE 



Develop analytical and ground test methods for predicting 
0-g performance. 

Develop techniques for system identification on orbit. 
Develop controllers using on-orbit system identification 
data. 
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HARDWARE DESCRIPTION 
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Making a precision spacecraft operate as required once on orbit implies 

TWO APPROACHES. FIRST, THE CONTROL CAN BE BASED UPON ANALYTICAL PREDICTIONS OF 
0-G BEHAVIOR. SECOND, THE BEHAVIOR OF THE SPACECRAFT CAN BE MEASURED ON ORBIT 
AND THE CONTROL CAN BE REDESIGNED. MACE FOLLOWS BOTH PATHS IN ORDER TO 
DETERMINE HOW MUCH PERFORMANCE CAN BE ACHIEVED THROUGH ANALYTICAL 
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TO DERIVE CONTROL FOR IMPLEMENTATION ON ORBIT. BY IMPLEMENTING THESE 0-G 
CONTROLLERS AT VARIOUS LEVELS OF CONTROL AUTHORITY, DEVIATIONS BETWEEN 
MEASURED AND PREDICTED PERFORMANCE REVEAL LIMITATIONS IN PREDICTION 
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SCIENCE APPROACH 
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SUSPENSION AND GRAVITY EFFECTS 
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Finite element models of both 0-g and i~g behavior have been compared 

WITH GROUND TEST MEASUREMENTS. THE 1-G FINITE ELEMENT MODEL PREDICTIONS OF 
TRANSFER FUNCTIONS WHILE INVOLVE MOTION OF THE TEST ARTICLE IN PRIMARILY THE 
VERTICAL PLANE COMPARE WELL WITH THE MEASURED TRANSFER FUNCTIONS. 




MACE 1-G AND 0-G MODELS 
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A CLOSED-LOOP ANALYSIS OR GRAVITY PERTURBATIONS WAS ALSO PERFORMED. A 
MODEL OF 0-G BEHAVIOR WAS USED TO DERIVE CONTROL. THE SPACECRAFT PERFORMANCE 
WAS THEN EVALUATED BY IMPINGING THIS 0-G CONTROL ON BOTH THE MODEL OF 0-G AND 
1-G BEHAVIOR. AS MIGHT BE EXPECTED, WHEN THE 0-G CONTROL IS EVALUATED ON THE 0-G 
MODEL, PERFORMANCE IMPROVES WITH INCREASING CONTROL AUTHORITY. HOWEVER A 
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GRAVITY INFLUENCE ON CONTROL 
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Performance vs. control authority Performance vs. control authority 
for entire model. for flexible portion of model. 


Measurement based control experiments have also been conducted on the 
ground, three different control topologies were used. 



Space Engineering Research Center 



. T3 

S § 
?— 
si 

<X> — E 
n a .s 
•g XI 5b 
6 g o 

§ ‘5b *3 

£ S« 

3 as 2 


be 

C a$ 

*0 | 
0) U* 
$ U 
Ch Q 

<U -9 

'bb > S 
g *Sb 

•a* 

• f* CQ 

U ^ 
c c 

• m a 

*■0 -2 
cd « 

,l/-e 

§, c 

I hN -FN 


S « 

ja u 

•i i 

be ^ 

c 3 

^ s 3 

<W g, 'w 

cS ? § 
S’ 

JU ja ^ 

'bb 6 13 
S’Ee -5 






Spate 'Engineering %{seartft Center 



83 a 

I ® 


5 * 
2 ° 


03 03 
03 ft 

^ s 

s* 

» a 

g e 

2 1 
~ 2 
H fc 
03 w 
5fa 55 
O O 

hH I— I 

fa H 

H § 

S3 s 
^ fa 

3 « 


Q 03 

< fa 
O fa 

j co 

£ 5 

< 3 


E o 
H o 

fa ^ 

O Q 

- , fa 
H co 
fa O 
O fa 

3 ° 

j E 

< ^ 

i— i ^ 

ft £ 
03 O 
g S3 
55 co 


£ fa fa 
O S3 S3 

O* P ^ 

i-5 03 q 
w E E- 

3 § O 

« ^ 55 

CO ^ 3 
CO < E 
S3 O z 

<J M O 
r K £ fa 
U Z CO 

° H 

2 Q 03 
g S 03 

« | ® 

•= 1 1 
<*3 S 


O < S3 

O « Eh 

h-3 Z fa 
A co O 

fa g fa 
co fa o 
O o Z 

o S 1 
u « § 
03 S3 « 
S3 en o 

03 ° fa 
H S fa 
ft S fa 

ip 

SIS 
2 £ * 


a « 

&H W 

z h 
o 3 
S 

fa ft 
Z 5 
o O 


5>h 03 
CQ fa 


a 3 

03 O 

£ 3 

fa M 

g a 
fa o 

CO 

H Q 
CO £] 

13 I 


W Cu 

si 

p 

g H 
S fa 
S3 o 
u 

_ w 

Si 

£ S3 

03 § 
O g 
03 D 
05 O 


3 £ 

O O 


O 'jj 

fa 

S S3 

Eh O 
t* H - 1 

3 S3 


Q e" 
< . 
O E-I 


< S3 

fa 03 


eg 

o fa 


w e 

u ^ 
Z fa 

fag 

03 S 

B s 
^ a 
s£ 


s a 

or « 

03 E-> 

0 w 

8 a 

fa a 

<3 2 

« 5 

1 s 

O 03 

fa £ 

22 3 
o w 
z 6 


§ ^ fa 

a fa S3 

•fa 5 ^ 

c z 

°|i 

2 5 & 

£j fa 

> 03 Q 

03 W £ 


r h H 

fa fa U 
fa fa 2 

fa Z ^ 

fa g S 
fa ^ s 
w w * 

fa CO co 
fa fa < 

Q g > 

2 fa fa 

t P s 

H a o 

03 § a 

o g w 

fe Ew 

co 5 6-1 
Z < fa 

O fa o 

h O rv 

ft g 


S3 O < 
fa O Q 


g S fa 

IP 

a fa s 

age 

a 3 g 

3<S 

£ O I 
fa ^ S 

gg§ 

fa Q w 

fa fa O 
fa fa Eh 

ill 

r_*i 'Z fa 


Eh W 

0 *fa 

1 3 

05 a 
fa a 

Q fa 

fa a 

% » 
fa o 

o 53 

E S 

z ® 
o o 
fa o* 
W t-S 
J J 

9 g 

s fa 
M O 

6 a 


^ Eh 

0 fa 

* 5 

CO fa 

a ^ 

S3 

j B 

1 < 

1 ^ 

8 B 

ss « 

is 

Oh Tf 

O co 

w w 

fa Q 
55 O 

< a 

S < 
O p5 

3e 


fa Eh W 

fa 55 ffi 

Q 03 £ 
O 05 p. 

s e s 


w a 

fa Q 


Oh CO 

a 03 


P ^ 2 

fa tf H 

fa fa fa 

s > <y 


fa fa fa 

s > ^ 

03 fa 
w fa « 
S3 03 
HOW 
fa S3 
fa W *- 
fa fa ffi 

9 9^ 
2 S g 

w w g 
a m w 

fa ^ 

a s a 

fa S3 S 

w S > 


E o 
H W 

< 03 

Q O 
Q fa 
O w 

S S 

Is 

fa M 

a a 


ft a 

fa Q 


Eh <3 


E5 p 3 

P|| 

g|B| 



Eh 

° S 

£ o 


5S”5 

BlSg 

W 05 O fa 

S3 03 £ r- , 

O h H W 
hh <; co a 
a fa fa P 


B|*l 

g g 3 1 

H a 03 w 
w o 9 fa 
fa w 2 2 



SINGLE-INPUT \ SINGLE-OUTPUT CONTROL 
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IN ORDER TO REALIZE THE SECOND AND THIRD CONTROL TOPOLOGIES, TWO ADDITIONAL 
TRANSFER FUNCTIONS NEEDED TO BE MEASURED. AGAIN, THESE TRANSFER FUNCTIONS 
WERE FITTED USING A NONLINEAR FREQUENCY DOMAIN FITTING ROUTINE WHICH 
MINIMIZES THE DIFFERENCE BETWEEN THE LOGARITHMIC VALUES OF THE COMPLEX 
TRANSFER FUNCTION DATA AND FIT. THIS ENABLES THE ROUTINE TO FIT THE POLES AND 
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SINGLE-INPUT \ TWO-OUTPUT MODEL 



Measurement model and transfer Measurement model and transfer 

function data from gimbal torque function data from gimbal torque 

to gimbal relative angle. to bus inertial angle. 




In the second control topology, the lack of phase margin near 9 Hertz 

GENERATED A NEED FOR LEAD COMPENSATION. A PAIR OF ZEROS BELOW 9 HERTZ FOLLOWED 
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SINGLE-INPUT , TWO-OUTPUT WITH 
PAYLOAD AND BUS PENALTY 
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FLIGHT EXPERIMENT RESOURCES 





Test article and support equipment stored in 3 middeck lockers. 
ESM stored in a fourth middeck locker. 

Wiring is pre-integrated in the test article for ease of assembly. 



Self-explanatory. 




i j 



Space Engineering Research Center 


FLIGHT EXPERIMENT SCIENCE OPERATIONS 
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performance can be tuned on orbit. 


Self-explanatory. 
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over an order of magnitude improvement in pointing accuracy. 





THE MIDDECK ACTIVE CONTROL EXPERIMENT: 
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GRAVITY AND SUSPENSION INFLUENCES 
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absence oi gravity and a suspension system. 

To develop rule-of-thumb predictions for the magnitude of the various 
gravity effects based on beam equivalent approximations of the 
suspended structure. 
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SUSPENSION EFFECTS BECOME 


APPROACH 
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become important for suspended structures. 
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OVERVIEW: Gravity and 
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Perturbations to the structure (i.e. the A matrix) can result in 
perturbations to the sensor and actuator performance 
depending on their location and the extent of the perturbation 
to the structure. 
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MODELLING OF GRAVITY AND SUSPENSION 

EFFECTS ON STRUCTURE 
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Suspension System 

The effects of the suspension system on the structural dynamics 
are captured by including the suspension system in the system 
model before incorporating the effects of gravity. 
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The effects of gravity on the B and C matrices have been identified and 
non-dimensional sensitivity measures show that these effects are 
especially important at low frequencies but are only significant with 
near-horizontal oriented devices and at those points where the 
rotations are large with respect to the displacements. 
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Similar expressions hold for the effective input of a PMA mounted to 
beam in bending and torsion vibration. 



APPLICATION TO MACE 
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MACE EM CONFIGURATION STUDY 
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EIz/EIy = 1/3 - rectangular section struts, destiffened about z. 
Out-of-plane - performance payload is swung 45° out-of-plane. 
Flex. App. - flap-type flexible appendages added to node 1. 
L-shaped - downward 90° bend is put in bus at node 2. 
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A gravity and suspension system modelling procedure has been 

developed which captures suspension effects, stiffness effects, static 
pre-deformation effects and direct sensor and actuator effects. 
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LQG controllers designed for on-orbit operation but tested on the ground- 
based model were shown capable of first destabilizing suspension 
modes and then flexible modes at higher levels of control authority. 
This makes ground testing of candidate on-orbit controllers difficult. 
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Identification For Robust Control 
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Three Levels of Identification 
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1. IN THE RICCATI EQUATION CONTROL PROBLEM THE OBJECT TO BE CONTROLLED IS 

A RICCATI EQUATION FOR THE AUGMENTED COVARIANCE MATRIX (INCLUDING BOTH 
STATES AND MODEL PARAMETERS). THE CONTROL IS INPUTS FOR IDENTIFICATION TAKING 
INTO ACCOUNT CONSTRAINT ON STRUCTURAL EXCITATION. THE CRITERION IS ROBUST 
CONTROL PERFORMANCE. 
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PETERSEN-HOLLOT'S BOUNDS ARE MODIFIED FOR RANDOM CORRELATED 
ITERS. THE CORRESPONDING MODIFIED RICCATI EQUATION IS DERIVED. 
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b). Duality principle for design of dynamical feedback 






• Realistic statistical model of uncertainty 
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Advantages of " Post-ID " Model of Uncertainty 
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PRACTICAL REALIZATION OF THE ABOVE APPROACH FOR MACE CONSISTS OF THREE 
STAGES: 1) SIMULATION 2) GROUND EXPERIMENT AND 3) EXPERIMENT IN SPACE. 
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ROBUST MULTIVARIABLE CONTROL WITH 
PARAMETRIC UNCERTAINTY 
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MULTIVARIABLE IDENTIFICATION 
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Methodology tested on interferometer testbed using accelerometer 
measurements. Extended to laser measurements (R. Jacques, F. 
Aquirre, and others). Will be used by Douglas for MACE. 
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onlinear method with new cost function (Jacques): 

J = 11 logG- log G m \\ 2 
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THESE PLOTS ARE THE SAME AS THE PREVIOUS ONES, EXCEPT THAT THE FREQUENCY IS NOW 
IN A LOG SCALE. THIS EMPHASIZES THE LOW-FREQUENCY PORTION OF THE GRAPH THE FIT 
AT LOW FREQUENCY IS NOT AS GOOD AS THE HIGHER FREQUENCIES BECAUSE THERE WERE 
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MACE EXAMPLE: 
TORQUE WHEEL TO ENCODER 
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MIMO STAtE SPACE MODEL 
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FUTURE WORK 
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January 23, 1992 



Sensor and actuator placement for robustness is the subject of a PhD thesis by Eric Anderson 
This should be completed by Feb. 1993 . 
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Outline 




Ml‘T Space ‘Engineering ‘Research Center 


The sensor/actuator selection problem has been recognized for a long time. Although there are 
many approaches available, most share a common attribute: the selection is based on a certain 
model. Notably, actuator and sensor noise have been considered (DeLorenzo, et al.), and costs 
based on the sensitivity of performance to uncertain parameters developed (Skelton). 
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Sensor! Actuator Selection 
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An explicit inclusion of the model uncertainty is included in robust control designs. It is 
believed that the uncertainty should also be represented in the sensor/actuator placement 
problem. 




Sensor! Actuator Selection for Robustness 
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The plot is a qualitative speculation of the possible performance improvement. The 

perfomance objective (y-axis) is to be minimized in the presence of some number of uncertain 
parameters (x-axis). 
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Achievable Performance Robustness 
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System Design Methodology 
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The description of the uncertainty is critical. The errors typical of finite element modeling of 
complex built up structures must be represented. Regardless of the later ability to do in-line 
identification of control loops, an inaccurate model will be used to place the transducers. 
Further, any ground test of a large structure (if this is even feasible) may not accurately 
capture 0-g dynamics. 
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Representing the Uncertain System 
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Refer to Anderson PhD for much more complete description of the costs. To date, these results 
are very preliminary. It is important however to recognize the need for closed loop evaluation. 
A topic of current interest is whether one criterion can be used for the s/a selection and an 
entirely different one for the control design. 
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Design Algorithms 
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Approach 

Structural dynamic analysis 
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distribution function, and average 
Check controllers on evaluation model 


The primary test case for the study is the interferometer multipoint alignment testbed. 
However, a much simpler beam problem is being considered first. In this case, the number of 
possible transducer combinations is easily within brute force direct search techniques. More 
sophisticated search algorithms, restrictions, and approximations will have to be applied to the 
placement problem 
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Simultaneous Sensing and Actuation 
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Physical Interpretation 
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This is the laser performance outputAlisturbance force input for the LQG controller. The gain 
was kept low because of low frequency unmodeled dynamics of the optical bench. 

Strain rate feedback (SRF), positive position feedback (PPF) and LQG compensators were all 
designed based on the system model and implemented on the beam system. Results agree well 
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Expression for envelope amplitude at tN, 
the time of the final impulse. 
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WE CAN EXTEND THIS METHOD TO 
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These four equations are repeated for each mode 
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AN EXAMPLE SOLUTION OF IMPULSES 
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THREE STEPS, AS SHOWN IN THE FIGURE. THE INPUT SHAPER CAN, THEREFORE, BE 
USED IN REAL-TIME, AS A FILTER TO DESIRED INPUTS. 
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DISCOS MODEL OF MACE. 



Body 1 = Bus + Gimbals + Torque Wheels 
Body 2 = Payload 1 (Rigid) 

Body 3 = Payload 2 (Rigid) 
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ENDPOINT RESPONSE TO UNSHAPED SLEW 
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PRESSURE ACTUATOR WITH 
VISCOUS FLUID DAMPING 
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anisotropic (composite) cavity to maximize axial stroke and 
hoop stress tolerance 

-- Incorporation of actuator into a controlled truss structure 
STATUS: Initial component design completed 
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Passive Control/Damping: Spangler, Hall 
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* Piezo is both sensor and actuator, truly collocated. 

* G22 is an electrical impedance which exhibits the com- 
bined effects of the piezo capacitance and the structure’s 
dynamics. 
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The approach to control design involves defining a cost functional of z 

SUCH AS ITS H2 OR Hoc NORM, AND OPTIMIZING IT OVER THE SET OF ALL PASSIVE 
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Passive Control: Spangler, Hall (cont.) 
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Develop a more fundamental, elegant solution technique 
for the W2/RC design problem (riccati equations). 

Experiment on the Interferometer Testbed. 
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Extensions to tracking control and to controllers with robustness to 
parameter uncertainty seem possible. 
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